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ABSRACT 


Predictions  of  the  sound  pressure  field  in  a  wedge-shaped  ocean  overlying  a  fast 
bottom  have  been  generated  by  computer-implemented  calculations  based  on  a  method 
of  images  model  and  two  parabolic  equation  models.  Comparisons  of  the  results  show 
significant  agreements  and  also  some  disagreements.  The  geometric  shape  of  the  wedge 
for  the  image  calculation  consists  of  a  plane,  sloping  penetrable  bottom  and  a  plane, 
horizontal  pressure  release  upper  surface.  This  shape  models  a  real  continental  shelf. 
The  shape  used  in  the  parabolic  equation  models  is  a  pressure  release  plane  surface 
above  a  conical  bottom.  Two  PT  models  were  studied.  One,  from  a  literature,  is  based 
on  a  wide-angle  parabolic  equation  arid  providing  a  contour  plot  ol  TL  as  a  function 
of  position.  The  other,  resident  on  the  NPS  computer,  is  based  on  an  implicit  finitc- 
differencc  (I I'D)  algorithm.  The  results  show  that  there  is  fair  agreement  among  the 
different  models.  Comparison  of  the  image  method  and  the  second  PL  model  shows  a 
2  to  3  dB  difference  in  transmission  loss  near  the  surface  and  good  agreement  deeper. 
The  transmission  loss  predicted  by  the  first  PL  model  differs  from  both  other  models  by 
4  to  5  dB  near  the  surface  and  at  the  middle  depths.  Near  the  bottom  they  all  agree  well. 
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I  INTRODUCTION 


Acoustic  propagation  in  an  ocean  with  a  sloping  bottom  is  of 
growing  concern  both  for  theoretical  reasons  and  because  of  the 
importance  of  knowing  the  performance  of  acoustic  sensors  located 
over  the  continental  slope.  This  problem  has  been  studied  by  many 
scientists.  In  this  research  we  compare  and  contrast  the  results  of 
three  different  models,  predicting  the  pressure  field  in  a  wedge- 
shaped  ocean  with  a  pressure  release  surface  and  an  acoustically  fast 
bottom.  One  image  method  and  two  PE  models  were  selected. 

The  computer  model  for  the  image  method  was  written  in  Basic 
for  use  on  a  PC  and  in  Fortran  for  use  on  the  NPS  IBM  3033  computer 
main  frame.  The  majority  of  the  computations  for  the  image  method 
was  done  using  the  Fortran  version. 

A  PE  model  contained  in  "Ocean  Exec”  library  on  the  NPS  IBM 
3033  computer  was  tested,  but  found  to  be  too  cumbersome  and  time- 
consuming.  Modeling  a  sloping  bottom  with  acceptable  accuracy 
required  100  pairs  of  range-depth  input  which  is  tedious.  In  addition 
the  program  needs  knowledge  of  certain  parameters  for  a  basement 
layer  which  is  required  to  exist  below  the  bottom.  The  output  gives 
transmission  loss  at  a  fixed  depth.  To  find  the  TL  as  a  function  of 
range  for  various  depths  between  the  bottom  and  the  surface 
necessitated  running  the  program  many  times  which  would  have 
required  an  unacceptable  amount  of  computer  time.  For  these  reasons, 
this  PE  model  was  not  used  for  this  work. 
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Two  other  PE  models  were  found  to  be  more  acceptable.  One  was  a 
model  used  at  SACLANT  by  Tindle  and  Jensen  (Ref.  1]  which  we  will 
refer  to  as  PESAC.  The  other  was  a  model  based  on  an  implicit  finite- 
difference  (IFD)  algorithm  implemented  at  NPS  by  Jaeger  [Ref.  2);  we 
refer  to  this  model  as  PEIFD  in  what  follows.  PESAC  was  not  available 
for  use,  but  some  of  its  output  was  in  the  referenced  article  in  the 
form  of  a  contour  plot  of  TL  vs  position  in  and  below  the  wedge  for  a 
stated  set  of  water  and  bottom  parameters.  The  PEIFD,  being  available 
at  NPS,  was  used  to  generate  output  for  the  same  set  of  parameters. 
The  results  of  the  image  method  model  and  the  two  PE  models  are 
compared  in  the  analysis  that  follows. 
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IL  DEVELOPMENT 


A  IMAGE  THEORY 

In  1978.  with  the  aid  of  a  computer  model  based  on  the  method  of 
images,  Coppens,  Sanders,  Ioannou,  and  Kawamura  (Ref.  3]  predicted 
and  measured  the  pressure  amplitude  and  phase  in  the  upslope 
direction  along  the  bottom  of  a  wedge-shaped  fluid  layer  overlying  a 
fast  fluid  bottom.  Baek  [Ref.  4]  used  the  same  model  in  1984  to  predict 
pressure  amplitude  and  phase  everywhere  within  the  wedge  in  the 
upslope  direction.  In  the  same  year,  LeSesne  [Ref.  5]  implemented  a 
n  odel  developed  by  Coppens  and  Sanders  which  is  not  limited  to  up 
or  downslope  direction.  In  the  method  of  images,  boundaries  are 
replaced  by  images  of  the  source,  placed  geometrically  to  represent 
the  propagation  path  between  source  and  field  point  and  given 
amplitudes  and  phases  to  reproduce  the  effects  of  interactions  with 
the  boundaries.  For  a  wedge-shaped  duct,  images  lie  on  a  circle 
centered  at  the  apex  of  the  wedge.  The  source  and  each  of  the  images 
radiate  spherical  waves  of  appropriate  phases  and  amplitudes.  The 
phase  coherent  summation  of  these  waves  yields  the  total  pressure 
and  phase  at  any  field  point  in  the  wedge.  In  the  computer  program 
the  number  of  images  is  made  finite  by  insisting  that  the  wedge  angle 
P  is>  an  integral  submultiple  of  180°.  Then,  the  number  of  images 
(including  the  source)  in  upper  or  lower  half-spaces  is 

N  =  (180/p)  (1) 
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Figure  1.  Image  structure  for  a  wedge  shaped  duct  (side  view) 


The  geometry  of  the  image  method  is  represented  in  Fig.  1.  Figure 


2  contains  the  definitions  and  symbols  for  source  and  receiver 
parameters. 


surface 


$  =  wedge  angle 

Ri  =  scaled  distance  from  apex  to  the  source 

R2  =  scaled  distance  from  apex  to  the  receiver 

Yo  =  scaled  distance  along  the  shore  between  the 

receiver  and  the  source 

y  =  source  angle  measured  upward  from  the  interface 

5  =  receiver  angle  measured  upward  from  the  interface 

Di  =  the  ratio  of  medium  density  to  the  bottom  density 

OC  =  the  ratio  of  medium  (water)  sound  speed  to  the 

bottom  sound  speed 

0C  =  critical  angle  for  the  bottom 

X  =  scaling  distance 

a/k2  =  loss  term  in  the  bottom  (used  in  the  reflection 
coefficient) 


R(n)  =  range  from  the  nth  image  to  the  receiver 
Figure  2.  Definitions  and  notations  of  source  and  receiver  parameters 
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For  the  wedge-shaped  duct,  sound  will  propagate  from  the  source 
to  the  receiver  through  paths  that  can  be  divided  into  two  families:  (1) 
direct,  surface  reflected,  surface  then  bottom  reflected,  surface- 
bottom-surface,  and  so  on;  (2)  bottom  reflected,  bottom-surface 
reflected,  bottom-surface-bottom,  and  so  on  For  example,  the  bottom- 
surface-bottom  path  to  the  receiver  in  Fig.  3(a)  can  be  interpreted  as 
the  straight  line  connecting  the  3rd  image  with  the  receiver  in  Fig. 
3(b). 

The  angle  0n  of  the  nth  image  measured  from  the  bottom  is  shown 
in  Fig.  1, 


=  (J(n-l)  +  y. 

for  n  odd 

(2.1) 

=  pn  -  y. 

for  n  even 

(2.2) 

The  pressure  from  an  image  can  be  determined  knowing  the 
distance  between  the  image  and  the  receiver  and  the  total  reflection 
coefficient.  The  total  reflection  coefficient  is  the  product  of  all 
reflection  coefficients  along  the  path  from  image  to  receiver.  With 
pressure  one  meter  from  the  source  normalized  to  unity,  the  pressure 
at  a  field  point  from  a  specific  image  is  found  by  dividing  this  total 
reflection  coefficient  by  the  distance. 
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source 


°  receiver 


APEX 


(b) 

Figure  3.  Visualization  of  (a)  real  path  and  (b)  equivalent  path  in 
image  method 


The  distance  between  the  receiver  and  the  nth  image,  as  shown  in 
Fig.  4  is 

R(n)  =  VRi  +  r2  ~  2R!R2  cos(9n  -  8)  +  Y§  (3  ^ 

for  the  upper  group  of  images,  and 

R'(n)  =  +  r2  -  2R!R2  cos(9n  +  5)  +  Y§  (3  2) 

for  the  lower  group  of  images. 


From  Fig.  2,  each  image  interacts  with  the  bottom  and  surface  a 
number  of  times:  so  that  the  pressure  at  the  receiver  from  a  specific 
image  is  proportional  to  the  product  of  the  relevant  reflection 
coefficients.  To  find  a  reflection  coefficient,  the  angle  of  incidence 
must  be  determined.  The  greater  the  index  n  of  the  image,  the  more 
reflections,  and  one  must  find  an  incident  angle  for  each  of  these 
reflections. 
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,  Rj  sin(0n  -  mP)  +  R2  sin(mp  -  8) 
sm(enm) - j^j 


Figure  5.  The  angle  of  incidence  0nm  of  the  upper  nlh  image  on  the 
mth  plane 
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for  the  lower  images.  [The  notation  is  0  for  the  upper  images  and  0'  for 
the  lower  images.] 

The  formula  used  in  the  program  for  calculating  the  reflection 
coefficient  from  a  lossy  bottom  is  derived  in  the  Appendix.  The  spatial 
phase  term  for  each  image  is  calculated  from  the  distance.  The 
pressure  from  an  image  is  determined  by  multiplying  the  spatial  phase 
factor  with  the  total  reflection  coefficient  and  dividing  it  by  the 
relevant  distance.  The  summation  of  the  complex  pressure  from  every 
image  is  the  complex  pressure  at  the  arbitrary  receiver  position. 

In  the  pressure  calculation  for  upper  images,  it  is  not  necessary  to 
calculate  reflection  coefficient  for  the  first  and  second  image.  These 
two  images  don't  interact  with  any  bottom  plane.  The  reflection 
coefficient  for  the  first  image  is  1  and  for  the  second  is  -1. 
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The  resulting  complex  pressure  from  the  upper  family  of  images  is 


pu  =  X  ^exp(-jkR(n)).(-l]'n'2'  fme™) 

n_lR(n)  m=1  (5.1) 

riR(0nm)  =  1.  for  n  =  1.2 

and  from  the  lower  images  is 

Pl=  I  -^-exp(-jkR'(n))  (-l)(n/21  flRO'nm) 

n=1Kln)  m=o  (5.2) 

In  these  expressions  n  indicates  each  image,  and  N  is  defined  in 
Eqn.  (1).  The  index  m  indicates  each  bottom  and  the  number  M  is  the 
total  number  of  the  bottom  planes  the  sound  meets  during  propagation 
to  the  receiver.  It  can  be  seen  that 

M  =  [(n-l)/2]  (6) 

The  total  complex  pressure  distribution  P  is  the  sum  of  the  Pu  and  Pj; 

(7) 


P=Pu+Pi 


a  PE  MODEL 

As  stated  above,  two  parabolic  equation  (PE)  models  were 
compared  with  the  image  method.  The  first  was  the  PESAC  model 
from  Jensen  and  Tindle  (1987)  which  is  based  on  the  wide-angle 
parabolic  equation.  The  pressure  field  in  the  model  has  cylindrical 
symmetry  about  a  vertical  axis  through  the  source  position.  The 
sloping  bottom  therefor  forms  a  conical  hill  for  downslope  propagation 
and  a  conical  bowl  for  upslope  propagation.  A  line  source  was  used  for 
the  majority  of  cases.  In  one  downslope  case  they  used  a  simple, 
standard  Gaussian  PE  source,  with  a  -3  dB  half-width  of  35°,  a  wedge 
angle  of  10°  and  a  frequency  of  5  Hz.  The  contour  graph  for  this  case 
will  be  compared  to  the  results  of  the  image  method  and  a  second  PE 
model. 

The  other  PE  model,  PEIFD,  was  drawn  from  the  thesis  of  Jaeger. 
This  uses  the  implicit  finite -difference  (IFD)  equation  solution  method 
to  solve  the  PE.  This  method  is  unconditionally  stable  and  has  the 
capability  to  incorporate  horizontal  or  sloping  interfaces.  The  Fortran 
program  in  this  thesis,  like  the  PESAC  model,  specifies  the  receiver 
position  with  range  and  depth  rather  than  range  and  angle.  The  input 
data  used  in  this  Fortran  program  were  the  same  as  the  input  data  for 
Jensen  and  Tindle. 
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EL  COMPUTATION  OF  IMAGE  METHOD 

For  the  environmental  parameters,  the  water  has  density  px  and 
sound  speed  Cj  and  the  bottom  has  density  P2,  sound  speed  C2.  and 
attenuation  a/k2. 

To  verify  spherical  spreading  close  to  the  point  source,  source  and 
receiver  were  positioned  as  in  Fig.  6.  With  the  receiver  in  the  same 
angular  plane  as  the  source,  the  receiver  range  from  the  source  was 
varied.  The  results,  Fig.  7,  verify  that  the  pressure  near  the  receiver 
has  a  1/r  dependence.  In  Fig.  7(a),  the  source  is  at  a  scaled  distance 
Ri  =  2  and  the  multiplication  of  the  range  and  the  pressure  is  close  to 
1,  as  it  should  be.  The  relation  between  scaled  and  real  distance  will 
be  illustrated  in  chapter  4  with  a  specific  case.  The  very  first  point  at 
unit  range  is  inaccurate  because  of  calculational  round  off. 

Figure  7(b),  where  the  source  is  at  the  scaled  distance  of  Ri  =  40, 
shows  the  same  dependence  within  calculational  inaccuracy.  In  this 
graph  there  were  no  results  at  receiver  ranges  of  1,  2,  3  m  because  of 
data  overflow  when  source  is  this  close  to  the  receiver.  This  overflow 
is  due  to  the  value  of  the  R(n).  From  Eqn.  (3.1),  R(l),  the  distance 
from  the  first  image  to  receiver  (direct  path),  is  very  close  to  zero. 
This  near  zero  condition  is  a  denominator  for  the  pressure  calculation 
of  Eqn.  (5),  thus  causing  an  overflow. 
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SOUI 


receiver 


Figure  6.  The  positioning  of  the  source  and  the  receiver 


The  next  three  figures  show  dB  loss  vs  position.  For  Fig.  8,  the 
source  is  close  to  the  apex  and  the  pressure  field  is  relatively  simple, 
corresponding  to  the  excitation  of  mostly  the  lowest  mode.  When  the 
source  range  is  relatively  small  this  mode  will  dominate.  The  contour 
plots  in  Fig.  9  and  10  show  that  placing  the  source  further  from  the 
apex  excites  the  second  mode.  Placing  the  source  even  further  from 
the  apex  will  result  in  higher  and  more  complex  modes  being  excited. 
All  the  distances  in  Fig.  8,  9  and  10  are  scaled  distances.  It  is  easy  to 
show  that  the  nondimensional  cutoff  distances  of  the  first  few  modes 
are  approximately  3,  5,  7,  11,  . . 


in 


0  2  4  6  8  10 


RANGE  FROM  SOURCE  (M) 

(a)  source  distance  Ri  =  2 


in 


0  2  4  6  8  10 


RANGE  FROM  SOURCE  (M) 

(b)  source  distance  Ri  =  40 

Figure  7.  The  range  R(l)  dependence  of  the  pressure  in  the  image 
method.  Frequency  5  Hz,  Dl  =  0.5.  CC  =  5/6,  (3  =  10°,  y  =  5  =  5° 
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319NV 

Figure  8.  Excitation  of  mode  1.  p  =  10°,  Ri  =2  and  the  rest  of 
parameters  same  as  in  Fig.  7 
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RANGE  FROM  SOURCE 


1  I 

8  9 

319NV 

9.  Mode  2  in  a  wedge,  Ri  =  < 
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RANGE  FROM  SOURCE 


TV.  DISCUSSION 


A.  IMAGE  METHOD  AND  PESAC 

The  contour  plot  from  Jensen  and  Tindle  (1987)  is  shown  in  Fig. 
12(a).  The  input  data  for  this  contour  are  given  by  Table  1.  The  input 
data  for  the  image  method  was  chosen  to  be  the  same. 

The  contour  lines  in  Fig  12(a).  which  indicate  dB  loss,  show  that 
propagation  is  almost  entirely  in  mode  1.  In  this  figure,  the  frequency 
was  5.0  Hz,  source  depth  (SD)  was  50  m,  water  depth  was  100  m  and 
the  maximum  range  was  10  km. 

To  compare  this  with  the  image  method,  geometric  parameters 
from  the  PE  model  were  translated  to  parameters  for  the  image 
method  using  Fig.  11.  Source  distance  Ri.  source  angle  y,  and 
horizontal  distance  x  for  image  method  can  be  found  from  the  wedge 
angle,  water  depth  and  source  depth  of  the  PE  model.  From  simple 
geometric  calculation,  x  =  567.13  m.  y  =  4.96°.  Ri  =  569.33  m. 


TABLE  1 .  PARAMETERS  FOR  MEDIA  USED  BY  JENSEN  AND  TINDLE  ( 1 987) 


densityf^rrr11 

sound  speed 

attenuation 

water 

pi  =  1000 

cj  =  1500  m/s 

bottom 

p2  =  2000 

c2  =  1800 

aX2  =  1  dB 
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X 


source  depth  50  m 
water  depth  100  m 


Figure  1 1 .  Conversion  of  source  distance  and  angle  for  image 
method  from  Jensen  and  Tindle  (1987) 

To  input  ranges  into  the  image  method  program,  they  should  be 
converted  to  a  scaled  range.  To  convert  Ri  and  R2  into  the  scaled 
range,  the  following  formula  is  applied. 


2  sin  9C  tan  (3 

where  X  is  the  scaling  range 

This  scaled  range  is  used  in  the  image  method  program.  The 
relation  between  scaled  range  and  real  range  is 

R(scaled  range)  =  R(real)  /  X 
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The  critical  angle  (0C)  is  calculated  to  be  33.56°.  With  this  critical 
angle,  a  wedge  angle  of  10°.  and  a  frequency  of  5.0  Hz,  X  =  769.48  m. 
This  gives  a  scaled  distance  for  Ri  of  0.74.  A  scaled  distance  of  R2  = 
13.2  corresponds  to  the  map  range  of  10.154  km.  (The  distance  R2  is 
measured  from  source  to  the  receiver  directly  in  the  image  method 
but  horizontally  in  the  PE  model.) 

The  attenuation  in  the  PE  model  of  a?i2  =  1  dB  is  equivalent  to 
a/k2  =  0.0183  for  the  image  method  since 

a/k2  =  a^2/(8.7  *  2  rc) 

The  transmission  loss  contours  from  the  image  method  are  shown 
in  Fig.  12(b).  These  two  contours  show  very  similar  shape.  Figure  13  is 
an  overplot  of  image -method  results  on  the  contour  plot  of  the  PE 
model.  Near  the  surface  and  in  the  mid-depth  level,  the  image  method 
has  3  to  5  dB  less  loss,  while  they  agree  well  near  the  bottom. 
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B.  IMAGE  METHOD  AND  PEIFD 

The  PEIFD  model  is  designed  to  generate  transmission  loss  both  in 
the  bottom  and  in  the  water  along  horizontal  range  from  the  source  at 
a  given  depth.  The  image  method  program  was  run  in  such  a  manner 
as  to  yield  comparisons  between  the  image  method  and  this  PE  model. 

The  input  file  used  for  this  purpose  is  shown  in  Table  2.  The 
number  in  column  1,  row  1  is  the  frequency  in  Hz  and  column  3,  row 
1  is  the  receiver  depth  in  m  which  was  varied  to  probe  the  entire 
sound  field,  both  in  the  water  and  in  the  bottom.  Rows  3  and  4  define 
the  sloping  bottom  and  the  definition  of  the  remaining  terms  can  be 
found  in  Ref.  2. 

In  the  image  method,  to  find  the  pressure  as  a  function  of 
horizontal  range  at  constant  depth,  the  receiver  angle  must  be 
changed,  as  the  range  is  changed,  as  shown  in  Fig.  14. 


TABLE  2.  INPUT  DATA  FILE  FOR  THE  PE 


5 

50 

500 

0 

4000 

10000 

0.0 

0.0 

50 

1000 

25 

0 

100 

10000 

1763 

-1 

-1 

1763 

1.0 

0.0 

0 

1500 

1000 

1500 

1763 

1500 

8000 

2.0 

1.0 

1800 

6000 


The  graphs  shown  in  Fig.  15  through  22  indicate  the  results  of  the 
PEIFD  model  with  a  solid  line  and  the  image  method  with  a  dashed 
line.  The  PESAC  model  is  annotated  with  the  'o’  connected  by  dashes 
for  selected  graphs.  The  graphs  are  dB  loss  versus  horizontal  range  in 
km  at  a  fixed  depth  and  all  the  source  and  receiver  parameters  are  the 
same  as  used  in  the  Jensen  and  Tindle  model.  The  difference  between 
the  image  method  and  PEIFD  is  approximately  2  to  3  dB  in  shallow 
water  and  the  two  lines  converge  as  depth  increases.  The  PEIFD 
results  in  the  bottom  also  have  a  similar  shape  to  the  PESAC  results. 
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(aa)  ssoi  NOiivDVdoad 

Figure  15.  dB  loss  .vs.  range,  depth  =  50  m 
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Figure  16.  dB  loss  .vs.  range,  depth  =  100  m 
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Figure  17.  dB  loss  .vs.  range,  depth  =  200  m 
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Figure  18.  dB  loss  .vs.  range,  depth  =  300  m 
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Figure  19.  dB  loss  .vs.  range,  depth  =  500  m 
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Figure  20.  dB  loss  .vs.  range,  depth  =  800  m 
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Figure  21.  dB  loss  .vs.  range,  depth  =  1000  m 
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Figure  22.  dB  loss  .vs.  range,  depth  =  1500  m 
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V.  CONCLUSIONS 

The  geometric  implementation  of  the  ocean  wedge  and  its 
boundaries  by  the  image  method  and  the  PE  models  are  different.  To 
validate  the  use  of  these  very  different  geometric  models  it  is 
important  to  know  how  well  they  agree  with  each  other  for  specific 
geometries.  This  thesis  was  an  attempt  to  compare  the  image  method 
and  two  PE  models  for  propagation  direction  down  slope  and  with 
only  one  mode  excited. 

Between  the  contour  plot  of  the  image  method  and  PESAC  model, 
there  is  a  3  to  5  dB  difference  at  mid-depths.  This  difference  may  be 
due  to  different  normalization  techniques  or  to  the  basic  difference  in 
geometry,  but  this  hypothesis  cannot  be  verified  because  the  PESAC 
program  was  not  available.  The  plot,  however,  shr.^s  a  very  similar 
shape  for  the  two  models. 

On  the  other  hand,  when  the  image  method  is  compared  with  the 
PEIFD  model,  the  results  show  good  agreement  overall,  with  slightly 
less  agreement  at  shallow  depths.  This  would  suggest  that  the 
different  geometry  may  not  lead  to  different  TL's  in  this  range  of 
comparisons.  Overall,  however,  since  the  image  method  compares 
quite  favorably  with  both  PE  models,  all  three  are  apparently  good 
tools  for  underwater  sound  prediction  for  downslope  propagation  with 
only  one  mode  propagating.  It  is  still  an  unresolved  question  how  the 
various  models  will  compare  for  upslope  propagation,  cross-slope 
propagation,  and  the  source  positioned  to  excite  more  than  one  mode. 
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APPENDIX 


The  reflection  coefficient  is  given  by 


— sin0-  —  sin0t 

V  -  Pi  °2 

— sin0  +  —  sin0t 
Pi  C2 

where  0  is  incident  angle 

0t  is  transmission  angle 

subscript  1  denotes  water  and  2  bottom 


and 

sin0t  =^/l-(c2  /  ex)2  cos2  0 

then 

— sin0t  =  ^(Cx  /  c2  )2  -  cos2  0 
c2 

With  lossy  bottom 


s± 

C2 


Cl 


c2(l-Ja/MSt(1  +  ja/k2) 


and 

(f1)2  -  (cj  /  c2 )2(1  +  Ja  /  k2  f  =  (ex  /  c2)2  +  2 j(cj  /  c2 )2 a  /  k2 
c2 

makes 

-r^sinOt  =  yjlci  / c2)2  -cos20  +  2j(c1  /c2)2a/k2 

c2 
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Let 


4^-sin0t  =  yja-jb  =  (a2  +  b2)1/,4e^2 


where 


_ib 

a  =  (Cl  / c2)2  -cos20  b  =  2(Cl  /c2)2a/k2  <P"  a 


With  trigonometric  identity 


JO/2  _ 


1 


l  +  Vl  +  (b/a)2 

2Vl  +  (b/a)2  ^ 


Vl  +  (b/a)2-l 


2^l  +  (b/aY 


Vb2  -2 


+  a^  +a 


\  2Vb2  +  a2 


-j- 


Vb2 '-2 


+  a  -a 


2^b2  +a2 


Utilizing  this 


V2— sin0t  =  y  Vb2  +  a2  +  a  - 
c2 


2  2 
^  +az  -a 


— sin0-~VVb2  +  a^  +  a  +  j— ^-VVb2  +a2  - 


V  = 


_  Pi 


V2 


V2 


— sin0  +  -^=VVb2  +a2  +a  -  j-L^Vb2  +a2  -a 

p,  V2  JV2V 

This  form  of  the  reflection  coefficient  is  used  to  find  pressure 
field  in  the  program  code  for  the  image  method. 
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